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We provide ready-to-use time-domain gravitational waveforms for spinning compact binaries with 
precession effects through 1.5PN order in amplitude and compute their mode decomposition using 
spin-weighted —2 spherical harmonics. In the presence of precession, the gravitational-wave modes 
. . . {£, m) contain harmonics originating from combinations of the orbital frequency and precession 

00 ' frequencies. We find that the gravitational radiation from binary systems with large mass asymmetry 

, and large inclination angle can be distributed among several modes. For example, during the last 

■ stages of inspiral, for some maximally spinning configurations, the amplitude of the (2, 0) and (2, 1) 

(N ■ modes can be comparable to the amplitude of the (2, 2) mode. If the mass ratio is not too extreme, 

the 1 = 3 and 1 = 4 modes are generally one or two orders of magnitude smaller than the £ = 2 
modes. Restricting ourselves to spinning, non-precessing compact binaries, we apply the stationary- 
phase approximation and derive the frequency -domain gravitational waveforms including spin-orbit 
and spin(l)-spin(2) effects through 1.5PN and 2PN order respectively in amplitude, and 2.5PN 
order in phase. Since spin effects in the amplitude through 2PN order affect only the first and 
second harmonics of the orbital phase, they do not extend the mass reach of gravitational-wave 
detectors. However, they can interfere with other harmonics and lower or raise the signal-to-noise 
O , ratio depending on the spin orientation. These ready-to-use waveforms could be employed in the 

data-analysis of the spinning, inspiraling binaries as well as in comparison studies at the interface 
5^ ' between analytical and numerical relativity. 
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I. INTRODUCTION 



> 

m 

ff^ ■ Coalescing compact binaries made of neutron stars (NS) and/or black holes (BH) can produce gravitational waves 
' (GW) strong enough to be detected by ground-based interferometers, such as LIGO [l|, Virgo and GEO 3, oper- 
. ating in the frequency range 10-10'^ Hz. Moreover, supermassive BH binaries could be observed at lower frequencies 

^— I ' 10~^-10~^ Hz and up to cosmological distances by the proposed laser space-based antenna LISA For detection 

OO , purposes, matched filtering is applied to noisy data in order to extract any signals that match members of the template 

O- bank gig. 

^ I Gravitational waves produced during the long inspiral phase can accurately be modeled by the post-Newtonian 
(FN) approximation to general relativity [1]. As the BHs approach each other and their velocities increase, the FN 
expansion is expected to become less and less reliable. Late in the evolution, non-perturbative information contained 
in numerical-relativity (NR) simulations and FN-resummed methods @, as well as perturbation theory need to be 
taken into account in building analytical templates for inspiral, merger, and ringdown. In this paper, we shall limit 
the discussion to the inspiral phase of coalescing BHs. 

In constructing templates for detecting inspiraling signals, it is recommended to account for all physical effects 
which contribute significantly to the gravitational waveform. Those produced b y th e spins of the binary constituents 
are among the most important ones, especially for asymmetric compact binaries jlOl|. such as NS-BH binaries and 
BH-BH binaries with component masses {1711,1112) G [5, 15]M0 x [1,5]Mq. For detecting such systems, one ma y b e 
able to employ phenomenological methods which capture the essential features of spinning, precessing waveforms [12l | . 
However, parameter extraction (isl . 0, [lH would warrant the inclusion of as much information about the spins of 
the binaries as possible, so that one should employ physical templates [l^ [13, [13 &t the highest FN order for this 
purpose. 
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For non-spinning compact binaries, the GW phase evolution has been computed thr oug h 3.5PN order [T^, [2O, 
HH, m, [2^ and the and hx polarizations are available through 3PN ord er [2I. [25l . l26l. [27l. [2^. For spinning, 
precessing binaries, the GW phase evolution is known through 2.5PN order [2^, l30| for spin-orbit couplings, and 
through 2PN order [3l| for spin-spin couplings (spin(l)-spin(l) and spin(2)-spin(2) contributions have been obtained 
in Refs. [s^. Issj). Spin-orbit and spin(l)-s pin( 2) effects in the /i+ and hx polarizations were computed through f .5PN 
and 2PN order, respectively, in Refs. [3J, l35[. ^ More recently, the spin(l)-spin(2) contributions at 3PN order in 
the conservative two-body dynamics were found employing either effective-field theory techniques [3^ . [ssl . [soj or the 
Hamiltonian formalism of Arnowitt, Dcscr and Misner |40l.l4ll.[4^. Now, spin(l)-spin(l) and spin(2)-spin(2) effects at 
3PN order in the conservative two-body dynamics are also available |43| . For including those higher-order spin effects 
in the GW phase evolution and polarizations, the results [sol. |40|. [4ll. l43j| need to be extended to the non-conservative 
dynamics, notably to the GW energy flux. 

The importance of using templates that have amplitude corrections beyond the leading PN order (henceforth referred 
to as Newtonian approximation ^) was emphasized by different authors in the context of ground-based [3 14^, [4^, l47j 
and space-based detectors [H, [43, [13, [EH 112, HI, Ull , both for detection and parameter estimation. So far, the effect 
of spins and precession on parameter estimation was studied in Refs. [Tsl. IsBI. IsgL IstI . but those studies were limited 
to non-spinning and Newtonian GW polarizations [13, [2^, [2^ . 

In this paper we provide ready-to-use /i+ and hx polarizations in time domain for spinning, processing binaries 
through 1.5PN order. The actual computation of the gravitational waveform hij through 1.5PN order was done by 
Kidder [s^l, as well as Will and Wiseman jsH], but the ready-to- use /i+ and hx polarizations at 1.5PN order were 
only written explicitly for strictly circular orbits for which spins are aligned with the orbital angular momentum. 
Recently, Ref. [53 has obtained the time-domain GW polarizations for generic orbits through 1.5PN order in the 
binary's comoving frame. The /i+ and hx polarizations derived in the present paper for spinning, precessing binaries 
through 1.5PN order reduces to that of Refs. [3^. [ssj in the aligned case except for a few typographical errors which 
we correct. 

In view of future studies at the interface between analytical and numerical relativity [s^, [HO, [HH, [H, [6^, [HI, [EE, [bgI . 
\QT\. [ssI. [bqI. Ftp! we decompose the time-domain /i+ and hx polarizations in spin- weighted —2 spherical harmonics and 
compute the modes, him, to 1.5PN order. We then consider spinning, non-precessing binaries for which we derive 
the Fourier domain representation of the generated gravitational waveform within the stationary phase approxima- 
tion (SPA). We provide a very compact way of writing the Fourier transforms of and hx which can readily be 
used for data analysis, for comparisons with numerical simulations, or for building analytical frequency-domain tem- 
plates including inspiral, merger and ringdown [tH. Fz^. The impact of spinning, precessing waveforms for parameter 
estimation will be investigated in a future paper. 

The remainder of the paper is organized in the following way. In Sec. |TT]we draw the source and detector frames, 
and introduce conventions and notations. In Sec. Illll we provide ready-to-use /i+ and hx polarizations in time domain 
for nearly circular orbits. The polarization modes him with respect to the spin-weighted —2 spherical harmonics are 
derived in Sec.[lVl The features of the modes when spins are present is then discussed in Sec. [V] Section IVll focuses 
on spinning, non-precessing binaries. Wc compute there the Fourier domain waveforms with spin effects through 2PN 
order in the amplitude and 2.5PN order in the phase before discussing the main features caused by higher harmonics. 
Finally, we summarize in Sec. IVIII our main conclusions. Appendices |X] and |B] present the GW polarizations and 
modes for precessing binaries on nearly circular orbits through 1.5PN order for generic inclination angles, whereas 
Appendix [C] shows the PN coefficients of the center-of-mass energy and radiative energy flux for non-precessing, 
spinning binaries. Appendix ID] gives explicitly the frequency domain amplitude coefficients with non-spin terms to 
2.5PN and spin terms to 2PN order. 

II. SOURCE FRAME, POLARIZATION AND PARAMETER CONVENTIONS 

To obtain the GW polarizations, it is useful to express the gravitational strain tensor, hij, in an appropriate source 
frame. Next, one specifies an orthonormal polarization triad composed of the direction of propagation N and two 



^ Note that spin-orbit effects through 2PN order in the h+ and hx polarizations were calculated in Ref. [36l |. However, Ref. [30| pointed 

out that a few multipole moments were computed incorrectly there. 
^ Note that the leading PN order in the polarization amplitude is proportional to 1/c"* when one turns the fundamental constants on. 

However, being the leading term in a PN expansion, it is has become common to call it Newtonian. 
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polarization vectors P and Q which arc used to construct the GW polarizations from the strain tensor [73 



h+ = -iP'P' -Q'Q')h,,, (2.1) 
hx = ^{P'Q' +Q'P')Kj. (2.2) 



The gravitational strain measured by a detector is then given by 

h{t)^F+h+{t)+F^h^{t), (2.3) 

where F+ and Fx are the antenna response functions that describe the detector's sensitivity to the two different 
polarizations. For laser interferometers with arms at a right angle, such as the LIGO and Virgo detectors, the 
antenna response functions for a GW coming from the sky location (9, (p) in the spherical coordinate grid built from 
the arm basis, with polarization angle "0, are [73| 

F+ = i(l + cos^0) cos20 cos2Vi-cos^sin20 sin2-!/i, (2.4) 

Fx = i(l + cos^^) cos 20 sin 2-0 + 008^ sin 20 008 2-0. (2.5) 

Note that the strain measured in a given instrument, h[t)^ is the same regardless of convention, whereas the wave 
polarizations depend on the choice of polarization vectors. Different choices of P and Q give different polarizations, 
but there is a compensating rotation of the polarization angle if) so that h{t) is unchanged. Here, we follow the 
convention of Refs. [1^, H^, in which 

P-^^' Q = NxP, (2.6) 
|N X Jol 

where Jq is the unit vector along the initial total angular momentum of the binary. In the absence of precession, 
the Newtonian orbital angular momentum La? = x v (with r, v, and being the binary separation vector, 
velocity, and reduced mass, respectively) is parallel to Jq. In this case, P coincides with the ascending node where 
the orbital separation vector crosses the plane of the sky from below. In the presence of precession, P is still defined 
as N X Jo/|N X JqI, but it is not in general the point where the orbital separation vector ascends through the plane 
of the sky. ^ 

For our source frame, we construct an adapted orthonormal basis (x,y,z) (see Fig. [IJ. We take the z-axis to be 
along Jq and the direction of GW propagation, N, to lie in the x-z plane, tilted by an angle 6 from the z-axis towards 
the X-axis. We describe the direction of the Newtonian orbital angular momentum with the spherical coordinate 
angles (t, a), where l denotes the angle between the orbital angular momentum and the z-axis while a is the angle 
between the a;-axis and the projection of the orbital angular momentum onto the x-y plane. For precessing binaries, 
as these angles vary in time, one must solve the precession equations to find their evolution. Notice that this source 
frame is the same as used in Ref. [sj], and depicted in Fig. 2 of that paper. 

We also find it useful to define basis vectors for the instantaneous orbital plane. These vectors have an implicit 
time dependence through the angles (i, a), and rotate about L at as it precesses. Here are their components in the 
(x, y, z) source basis: 

Jo X Ln . . „ \ fn '7\ 

XL = — = (— sma , cos a , ) , (^-'j 

|Jo X Ln| 

Yl = Ln X xl = (^ cosi cos a , — cost sin a , sim ) . (2-8) 

As an initial condition, we take the orbital separation vector n to lie along x^, at initial time, i.e., h{t = 0) = XL(t = 0). 
Then, we define the phase <!>(<;) to be the cumulative angle between XL{t) and n(t). 

h{t) = XL(t) cos$(t)-hyL(i) sin$(i), (2.9) 



^ Note that Ref. [33 | chooses polarization vectors rotated by 7r/2 relative to ours. This results in an overall sign difference from our 
polarizations, as can be seen by making the substitutions P — > Q and Q — ♦ — P, or by noting that GWs are spin-2 objects and flip sign 
under a 7r/2 rotation. As mentioned, the polarization angle of this convention is then rotated by 7r/2 relative to ours. This flips the sign 
of the antenna response functions as well, and so the same strain 1 12.31 1 is measured by either convention. 
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FIG. 1: We show (i) our source frame defined by the orthonormal basis (x, y,z), (ii) the instantaneous orbital plane which 
is described by the orthonormal basis (xliYLiLn), (iii) the polarization triad (N,P,Q), and (iv) the direction of the total 
angular momentum at initial time Jq. Dashed lines show projections into the x~y plane. 



\{t) = -xt(t) sin$(t) +yL(t) cos$(t). (2.10) 

We thus see that the phase <J>(t) measures how n has rotated relative to the vector x^. However, for a precessing 
binary, x^ is itself rotating about Ln- This means that the total rotation of n about Ln can be decomposed as a 
rotation of fi in the comoving basis parametrized by <i>(t) times a rotation parametrized by a precession phase due to 
the movement of the orbital plane itself. In the non-precessing case we have Jo || Ln and $(t) is expected to become 
the standard orbital phase whose time derivative is the orbital frequency. However, when Jq || Ln holds, we cannot 
define $ = for n = Jq x Ln/|Jo x Ln|, and we set $ = at the ascending node N x Ln = P, where the orbital 
separation crosses the plane of the sky from below. Now, $ = at the ascending node is achieved for x^ = n and 
P = fi, hence a = tt, so that the non-precessing regime is reached in the limit where /, = and a = tt for all time. 
This is applicable when the spins of the bodies are aligned or anti-aligned with the orbital angular momentum (or in 
the non-spinning limit). The waveforms are then greatly simplified. 
We define the following mass parameters 



M = 


mi +1712, 


(2.11) 


V — 


mi 7712 


(2.12) 


(mi 4- 7712)^ ' 


M = 




(2.13) 


5 = 


mi — m2 
mi -1- m2 ' 


(2.14) 


V = 


\(^-'')- 


(2.15) 



They are the total mass M, the symmetric mass ratio v, the fractional mass difference S, the chirp mass Ai. The 
symmetric mass ratio is bounded according to < < 1/4 and the fractional mass difference satisfies —1<6<1. 

The spin of a rotating compact body is of the order S ^ ml Uspm with I ^ Gm/c^. If the compact body is maximally 
rotating, then Ugpin ^ c and S ^ xGm'^/c. In words, from the PN point of view, the spin is formally of order 0.5PN. 
By contrast, if the compact body is slowly rotating, then I'spin c, and the spin is formally of higher PN order, 
S ^ xGm^ Wgpin/c^ ^ 1/c^. Throughout the paper, we use geometrical units where G = c =1. Henceforth, we shall 
work with the spin vectors normalized by the component masses as 

Xn = -%, «=1,2, (2.16) 
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so that Ixnl < 1 for objects that obey the Kerr bound on rotational angular momentum. Wc also define symmetric 
and anti-symmetric spin combinations as in Ref. [sH, 

Xs = ^{Xi+X2) , (2.17) 

Xa = liXi-X2)- (2.18) 

III. READY-TO-USE GRAVITATIONAL- WAVE POLARIZATIONS FOR PRECESSING BINARIES ON 
CIRCULAR ORBITS THROUGH 1.5PN ORDER: SMALL INCLINATION ANGLES 

The expression of the strain tensor hij for generic orbits through 1.5PN order was derived in Refs. [l^, [sHi and is 
given by Eq. (6.11) of Ref. [ssj. In this section we compute ready-to- use polarizations in time domain through 1.5PN 
order within the adiabatic regime where the binary inspiral is modeled as a quasi-stationary sequence of orbits assumed 
to be nearly circular. By nearly circular, we essentially mean an orbit that would be exactly circular, with separation 
vector ro of constant radius rg, in the absence of spins and gravitational radiation. The perturbation of the separation 
dr of such a motion is assumed to remain small with respect to rg on timescales on which the radiation-reaction effects 
can be neglected. This can only happen when the precession angles are at most of the same order of magnitude as 
the relative corrections induced by the spins in the dynamical quantities. Now, the evolution of Sr is governed by the 
radial part of the 1.5PN perturbation of the force per mass unit given in Eq. (2.1) of Ref. [IJ. It turns out that this 
perturbation depends on the spin exclusively through the two projections (S„ • Ln) with n = 1,2 which are almost 
constant apart from remainders that will contribute at higher orders in our weak precession hypothesis. In order to 
write the equation for Sr, we project the relative acceleration a in the basis {n, A, ^n}- For the sake of convenience, 
we introduce an "orbital" -like frequency Worb, defined as tJorb ~ (v ■ A)/r. The closure relation yields the following 
decomposition for v and a: 

V = rn + uJoibrX , 

a = (f - rtj^rb)n + (rtOoi-h + 2fcjorb)A - rworb • Ln , (3.1) 

with f = dr/dt. Splitting a into an unperturbed part slq plus a perturbation Ssl and using the equations of motion, 
we find finally that Sr satisfies the equation Sr + oj'^Sr ~ const, where loq is the constant angular frequency of the 
background motion. A particular solution is given by a constant perturbation, Sr ~ const., whereas the homogeneous 
solution satisfies an harmonic oscillator equation independent of the spin. 

By making the particular choice of a zero homogeneous solution, we can always eliminate the oscillations of r 
that are not directly linked to the non-zero spins of the BHs. Based on these observations, we shall define precisely 
a nearly circular motion to be a perturbed circular motion whose homogeneous radial perturbation solution ((Sr)hom 
is zero, as it would be for an exactly circular motion. ^ Assuming such a dynamics for our binary system implies 
that both 6r and r = r^ + Sr must be constant, provided we neglect higher order spin terms and radiative effects. 
We can generalize nearly circular motions to the case where spin precession angles are arbitrary in the absence of 
spin-spin interactions. This is achieved by introducing the concept of spherical motion defined as a motion having a 
constant separation r. It immediately follows from Eq. p.l[) that the full (conservative) acceleration is still of the form 
—ojl^^^rn. Moreover, when radiation-reaction effects are neglected, the orbital frequency computed from the 1.5PN 
equations of motion keeps being almost constant [sil. [s^. even for precession angles that are no longer small. This 
can be seen [29j by noticing that the only possible non-constant terms in Worb at the 1.5PN order come from the 
leading spin contribution of the equations of motion, and thus, are of the form (S„ • Ln). Their time derivative reads 
(dLN/'^i'S„) + (LN -dSn/dt). The first term is zero due to the precession equation dSn/dt ~ r2„ x S„, while the second 
term is a higher order correction quadratic in spins because of the approximate conservation of Lat. The treatment 
of the spin-spin dynamics is more delicate. A possible way to proceed consists in averaging the time dependent spin 
contributions in o^orb over one orbital period [s^ . . 

Introducing the invariant velocity, 

V = (Aft^orb)'/' , (3.2) 



* Though this type of motion can exist and is more general than the spin-aligned or anti-aligned case, it docs not necessarily represent 
yet the most likely evolution to be observed. 
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we reduce Eq. (6.11) of Ref. [3a | to nearly circular orbits and expand it in powers of v with the help of the relativistic 
extension of Kepler's law linking Uoi-h and M/r provided by Eq. (7.1) of Ref. [sH]. Schematically, we obtain 



Dr. 



Q,, + Pi/'Q,, V + (piQ,^. + piQ|0) + (P^/'Q,, + p3/2Qtaii ^ p3/2Q 



TT 



where SO indicates the spin-orbit terms; the tail integral P'^^^Q'^'' given by Eq. (6.11c) of Ref. [1^ reads 



(3.3) 



)3/2Qtail 



7r(A'AJ' - Wh^) + 121n ( — ) A^^ri^) 



TT 



(3.4) 



vq being an arbitrary numerical constant reflecting the freedom in the choice of the radiative time origin. The symbol 
TT on the square bracket indicates the transverse trace-free projection in the plane orthogonal to the direction N 
of the observer. We remind the reader that the non-spinning contributions to Eq. (|3.3p are known through 3PN 
order 



Apart from the spins, there are four vectors that appear in the expressions for the P"Q,y's in Eq. (|3.3p . In the 
source frame constructed in Sec. Illli they have the following (x,y, z) components 

n = (— sin a cos $ — cos b cos a sin <I> , cos a cos $ — cos l sin a sin <I> , sin t sin $) , (3.5) 



(sin a sin $ — cos l cos a cos $ , — cos a sin $ — cos l sin a cos $ , sin l cos $) 



(3.6) 



N = (sin 61, 0, cos 6*) , 



(3.7) 



Ln = (sin I cos a, sini sin a, cost) , (3-8) 

where $ is the phase defined in Eq. (j2.9p that measures how n has rotated relative to the vector xl. As xl is itself 
rotating about Ln for a precessing binary, the orbital frequency, or the total angular velocity of n about Ln, is the 
angular velocity of the motion of the binary within its instantaneous orbital plane, plus a precession velocity due to 
the movement of the orbital plane itself. To derive the relationship between the phase <I>(t) and the orbital phase (or 
carrier phase), we compute the derivative of n{t) by means of Eqs. (j2.7p . (j2.9p . obtaining 



dn f da\ ~ f di 

— = — + cos t — A + — sin $ 

dt \ dt dt \ dt 



sm L cos 4> — — Ln ■ 
dt ' 



(3.9) 



By imposing Ln = n x v/|n x v| = n x dh/dt/\n x dn/dt\, and using Eq. p.Sp as well as Eq. (|3.8p . we find that the 
term proportional to Ln in Eq. (|3.9p must be zero. Thus, we have dh/dt = fr/r^ + r/r = (v • A)A, where v • A is the 
orbital frequency cjorb defined before Eq. (|3.ip . which may be now interpreted as the angular velocity with which fi 
rotates about Ln- Identification with Eq. (|3.9p leads to the relation 



the phase <&(<:) being simply the integral 



(3.10) 



4>(t) = / [uorhit') - cos L{t')a{t')] dt' . 
Jo 



(3.11) 



Due to the freedom in the choice of the time origin by the radiative observer, /ly depends on an undetermined time 
scale or, cquivalently, on an arbitrary reference orbital frequency loq. The constant luq is actually associated to the 
presence of gravitational- wave tails and appears solely in logarithms of the form ln(worb/wo). Such contributions may 
be absorbed in the orbital phase by a redefinition of $ into a shifted phase "if Through 1.5PN order in the shift, 
we can pose = $ — 2v'^ ln(u;orb/wo). By plugging Eqs. (|3.5p - (|3.8p into Eq. (|3.3p . taking the combinations given in 
Eq. (j2.ip . and collecting terms by powers in v, we obtain the waveform polarizations 



Dl 



H 



(0) 



H 



(1/2) 



H 



(1/2, SO) 



H 



(1) 



H 



(l,SO) 
+ ,x 



H 



(3/2) 



H 



(3/2, SO)" 



(3.12) 
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The Newtonian. 0.5PN and IPN order terms were already computed explicitly in Refs. [3J, |35| [see in particular Eqs. 
(B2), (B3) of Ref. HJ]], but as a series expansion of M/r rather than v. Let us list for the reader convenience a few 
typographical errors we found there. In Eq. (4.9d) of Ref. [s^l, the factor of (1/6) (149 — Qv) has to be replaced with 
(1/6)(149 — 36z/); in Eq. (B2c) must be changed to —Q+; in Eq. (B3c) the right parenthesis is missing in the 
expression (cos^ i s\v? a + cos^ a); at last, in Eq. (B3j) cd should be read as —cd. In Ref. [1^, Eq. (F14b) must be 
multiplied by 2>v] in Eq. (F20) —iv has to be replaced with +'iv\ there should be an overall minus sign in front of Eq. 
(F25c). If we re-expand Kidder's polarizations in w, and correct all the previous typos, we obtain complete agreement 
with both results through IPN order. ^ 

The lengthy expression for the GW polarizations can be reduced to a much more compact form by noticing (see 
also Sec. IVD in Ref. [13]) that in the limit S <^ L the angle t can be considered a 0.5PN order correction. This can 
be seen from 



|Jo X L| 

JoL 



(3.13) 



if we neglect radiation reaction effects, i.e., we assume Jo = J, and use J = L 
then replace sint and cost in /i+^x with their Taylor series expansions in t. 



Si + S2 and S„ = 0{l/c). We may 



suit 



COSi = 1 - - +0(6*). 



(3.14) 
(3.15) 



However, the assumption S ^ L becomes less and less reliable for smaller mass ratio binaries. In fact, as a first 
approximation, we have Sn/L = [ran / M)"^ Xn v / v with v = (GMworb/c^)^^'^- Thus, even if Sn ^ 0{\/c), L can 
become comparable to Sn when v is sufficiently small. Moreover, we have assumed Jq = J in Eq. (|3.13p . but the 
latter is not exact when radiation reaction is included, and it can be strongly violated in presence of transitional 
precession [TSj . For these reasons, though we have decided to list in this section the GW polarizations expanded in 
t, we display in Appendix |^ the full expressions for generic inclination angles. For the (.-expanded polarizations, we 
find 



Hf =-(ce^ + l)cos2(a + «'). 



H 



(1/2) 



V 5 So 



H 



+ I) + ^') - ^ ( C0^ -I- 1 ) cos 3(a + *) 



ce^ 3ce^ 19 / , llce^ 19\ , 
^ + ^ + - + ice' + - _ ) ^ ) cos 2(a + *) 



(3.16a) 
(3.16b) 



-[l-ce^][iv -1] cos4(a + *) 



(3.16c) 



C0 



5ce2 19 



ce^ ce" , 49 



— + ( — + — li^l cos(a + *) - 27r(ce^ + 1) cos 2(a -f ^f) 

192 16 64 V 96 8 96/ ' \ ^ > v» ^ ; v -r ; 



+ (5se 



8W 45ce2 657\ /81c/ 9ce2 225^ 

+ \v C0s3(q; + W) 

128 16 128 / V 64 8 64 ' ' v -r ; 



(1/2, SO) 



^se^(^l-c/^ (2v-\ \ cos5(a + vI/) 
: — 2t cqsb cos(a + 2^*) , 



(3.16d) 
(3.16e) 



^ It is also worth noting that Ref. |34|| sets the origin of phase to bo at a point referred to as the ascending node and defined to be the 
point where the orbital separation crosses the x—y plane. This is in fact the same as our phase origin, = Jo X Ln/|Jo x LnI, but 
to reduce the possibility of confusion, we do not call this point the ascending node. We reserve this term to mean the point where the 
separation vector crosses the plane of the sky from below. 
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ce[Xa+ Sx"s ]-se[xl + ^xl cos(a + ^)-cg[xl + 5x^ sin(a + 4-) 
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■ v'^L Sg 



cg[ -ce'^+4+[ 3cg'^ + -\ iy\ cos(a + 2^-) + ( - ^ - 1 + ( cg^ + 3 ) ) cos(3a + 21') 
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fllcg^ 23 



+ ce ( ^ - 16j^ ) cos(3a + 4*) - (x^ + Sxl) sin 1- 



27/ 
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16 16 
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+ ] cos(a + + ^ (^ce^ + 1 j cos3(a + 1-) - -^(ce^ + 1) cos(3a + 1-) 
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cos(a + 3*) 



+ i CgSg 
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- cos(a - 2*) + - cos(q! + 2*) 

2 6 



H 



(0) 
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(1/2) 



Hi'^ = v'cg 



u(3/2) 3 



-2ce sin2(Q; + ^) , 
vdcgsg 

W . 17 



9 3 
- sin 3(q; + *) + - sin(Q; + 



13 



5 Sg 



135C02 603 
64 ^ "64" 



g + 4c02 jj/ j sin2(Q; + 1*) + se^( ^ | + ) sin4(a + 1*) 



_ 5£e_ ^ ^ \ I gjj^/ + I*) - 47r sin 2(a + *) 
48 48 ' ' 

171 libcg^ 



(3.16g) 

(3.17a) 
(3.17b) 

(3.17c) 



32 32 



v I sin3(a + ^) 



+ 5S9 {^{2y-l)sg^]&mb{a + ^) 
^a/2,SO) ^ sin(a + 2^) , 



(X^ + Sx"s) cos(a + 1-) + ce ce(x: + -^X?) - (Xa + ^Xs) sin(a + ^) 



+ LV S 



^(3/2,SO)^^3 



9 9c0^ 



sin(2a + 34') 



se'sinl-- (^^ + i^sin(2a + «')+ ^ 4. 2 
(^2(5x^ + (2 - i^) x^) (1 + cos 2(a + vj/) 
+ ( ^ 5xa + - (^ + 4c2)i. ) x^ - I 2<5x: + ( 2 + (3 + icg')u ) xf) ) sin2(« + v]/) 



(3.17d) 
(3.17c) 

'cesin2(a + 4'), (3.171) 
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+ LV Sg 



ce ix: + SXs) - Mxl + Sxt) sin* + - 3cg' + 6 + - - + 9cg' sin(a + 2*) 



16 



9 



ce 



3 



3c0^ + sin(3a + 2*) + - - + 8ce 



2W 



sin(3a + 4^') 



L^v5 CgSg 



3 IT 2T 1 45 

- sin(a - sin(a + H sin3(a + *) sin(3a + *) H sin(a + 35') 

8 16 16 16 16 



-sg sin(a-2*) 



-sg siii(a + 2*) 
o 



(3.17g) 



where sg and cg are shorthand notations for sin 61 and cos 9 respectively. In Sec. IIVI (see Fig. [2]), we shall discuss 
typical variations of the inclination angle l depending on spin orientations and binary mass ratios. Note that whereas 
the terms of h';^'^'^^ linear in x„ depend on the first harmonic of the orbital frequency, those of depend 

(2 SS) 

on its zeroth and second harmonic, and so do the terms of quadratic in the spin components, although we do 

not use them here. We include these 2PN SS polarization corrections when constructing frequency-domain waveforms 
for binaries having their spins aligned or anti-aligned with the orbital angular momentum in Sec. IVII The harmonic 
dependence of the polarization corrections produced by the spins can be understood from the explicit expression for 
/ly shown in Eqs. (4.9c) and (4.9d) of Ref. [H or Eqs. (F15a)-(F15c) of Ref. [35)]. The IPN SO contributions are 
proportional to the components of the orbital separation vector, n, which are themselves proportional to sin<l> and 
cos so that H^'^'^^ depend on the first harmonic of the orbital phase. Next, the 1.5PN SO and 2PN SS contributions 

are proportional to products of the orbital separation or instantaneous velocity unit vectors, n or A, and to products 
of sin$ or cos$. These can be re-expressed in terms of sin 2$, cos 2$ or constant quantities independent of so 
that '^^-^ and H^'^^"^ depend on the zeroth and second harmonics of the orbital phase. Because the expressions 

for hij in Refs. (33.[35l| are expanded in (M/r), while we use an expansion in v = (Mworb)^^'^, one has to convert from 
one expansion to the other by using Eqs. (7.1) and (F20) of Ref. [3^. In doing so, the w-expansion gains additional 
1.5PN SO and 2PN SS corrections proportional to the Newtonian order term depending on the second harmonic of 
the orbital phase. The IPN SO term is left unchanged. 

Ready-to-use time-domain GW polarizations for spinning, precessing binaries through 1.5PN order in phase and 
amplitude can be obtained by solving numerically the following equations: (i) The spin precession equations [29l. [30j 



~dt 
~dt 



— -!- = Oi X Si 



^ = X S2 



where at 1.5PN order 



5/3 



3 v 3 ^ , ^ 
- + -T4<5) Ln: 



(ii) The evolution equation for the Newtonian angular momentum 

Ln = — (Si + S2) 

(iii) The equation for the orbital frequency 



'"orb 



96 



1 



743 
336 



19 



113 
12 



113 



■ 4mr 



(3.18a) 
(3.18b) 

(3.19) 
(3.20) 

(3.21) 



Integrating Worb yields the orbital phase $oib(i) = Jq oih {f ) dt' . The GW polarizations (|3.12[) through 1.5PN order in 
phase and amplitude are computed by solving numerically Eqs. (|3.1ip . (|3.18|) . (|3.20p . and (|3.21[1 . In order to compute 
the GW polarizations (|3.12[) through 1.5PN order in amplitude, but at the highest available PN order in phase, one 
should replace Eqs. (jXTg)) . (jX^ with Eq. (7.5) in Ref. js^ and Eq. (32) in Ref. [H, respectively. 



IV. GRAVITATIONAL-WAVE MODES FOR PRECESSING BINARIES ON NEARLY CIRCULAR 
ORBITS THROUGH 1.5PN ORDER: SMALL INCLINATION ANGLES 

Due to the spin-2 nature of GWs, it is convenient to decompose the waveform components in the dyad {(P + 
iQ)/'\/2, (P — iQ)/V2} with respect to an orthonormal basis of spin ±2 functions that are defined on the 2-sphcrc 
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and belong to an irreducible representation of SO (3). Most commonly, the complex polarization 



^ ^ - iQ^ - iQ'' ^ ,,,, 
h = hjk ^ = h+-ihy,, (4.1) 

is expanded into the set of spin-weighted —2 spherical harmonics. Like the standard spherical harmonics, these 
functions of the two angles of spherical coordinates are labeled by a pair of integers, say {i/m), with i > 2 and 
m < 1^1 . The spin- weighted —s spherical harmonics associated to any such pair are given by ^ [7y] 



\e,4>) = (-l)V^<™(^)e""^ (4-2) 



with 



fc=max(0, m-s) ^ i^V' /V J \ / \ / 

and the orthogonality relation holds 



/ 



dn -sY^^^io, (j)) _sr^'"'*(6i, 0) = , (4.4) 



where dfi = sm6d6d(f> denotes the element of solid angle and S^^ is the Kronecker symbol. The integration is 
performed over the unit sphere, so that < 9 < tt and < ^ < 27r. The mode expansion of the complex polarization 
(14.11) is then of the form 



+00 £ 

h{eA) = Yl E /i^m-2F'™(^,<^). (4.5) 
The GW modes him are extracted using the orthogonality property (|4.4p by means of the surface integral 

h,n = Jdn h{e, 0) _2i^'™* {e, <p) , (4.6) 

where the star on the spin-weighted —2 harmonic indicates the complex conjugation. Therefore, the calculation of him 
requires the knowledge of the polarizations /i+ and h^ for an arbitrary value of the azimuthal angle of the direction 
N. In Sec, mil we have computed /i+ and h^ only for cf) = 0; however, a specific choice of the x-axis orientation cannot 
be responsible for any information loss. Thus, we must be able to recover h(6, (j)) from the expression of h{9^ 0) alone. 

The quantity /i at a given point depends on a number of parameters, such as l or the spin variables, and can actually 
be regarded as a function of 9, (f> as well as a function of the whole set of parameters that possess a geometrical 
character. More precisely, we may write h{9,(j)) = /i(6', i, a, $, x^, Xnj Xn)- Let us now introduce the projection 
basis (x',y',z' = z) obtained by applying a rotation of angle cj) about the z-axis on the vectors of the original basis 
(x,y,z). Let us also associate to each variable of a primed counterpart, which is defined in the same way as the 
unprimed variable but refers to the new basis rather than the original one. For instance, a' denotes the azimuthal 
angle of the orbital angular momentum measured from the fixed vector x' instead of x. In particular, we have 6' = 6, 
(/>' = 0, l' = L, a' = a — (j). The phase defined as the angle (Ln x z,n) = (Ln x z',n), is not affected by the 
transformation: $' = The x' and y' spin components can be obtained from the 2-dimensional formula for a passive 
rotation of angle (j), that is 

X'l = Xl cos cj) + xl sin (j), (4.7) 
X'n = "Xn sin(/)H-x^ cos</>, (4.8) 

while the third component is left unchanged. With our conventions, the polarization vectors in the new basis remain 
equal to P and Q respectively. Therefore, the complex polarization is identical to that of the old frame. Moreover, 



Our definition of —sY^"^{6, 0) differs from that of Ref. [T^ by a factor (—1)'" so that oY^"^{9, (f>) coincides with the most broadly used 
definition of Y^™'(d,<f>); for the d-matrix d^,^_^_^{9), wc adopt the same convention as Landau-Lifchitz |77| . 
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by construction of the primed variables, the functional dependence of h is the same as before, meaning that h = 
h{e,(l},L,a,^,Xn,xl,Xn) = Ke' ,l' ,a' ,^l_,x'l^x'n)- Hence the important relation 

h = h{9, (/), i, a, $, xl,X^, Xn) = HS, 0, t, a - $, cos (jiXn + sin 0xL - sin^x^^ + cos (j)xl,Xu) , (4-9) 

where the function h{9, 0, l, a, Xn^ Xn: Xn) is given by Eqs. (|3.16|) . (|3.17p for the t-cxpandcd expressions or by Eqs. (jA2[) . 
for the fuU ones. 

At 1.5PN order, the GW polarizations decompose into a sum of 3 terms, hQ{6, l, a, $) + X]n=i 2 Xn-hn(6', t, a, $), 
which shows that /i may be written as 

hoie, « - 0, $) + ^ W'^Qknie, a - 0, $) + e-'^Cnfc:(^, a - ^, $) + x^LKi^, ^ « - 0, <&)] , (4.10) 

n=l,2 

with ^„ = (x^ + ix^)/V2 and fc„ = (h^ + ihf^)/^/2. Each mode hem splits accordingly into 7 contributions: the spin- 
free term and 6 terms proportional to each of the spin variable components. These contributions arc parametrized 
by a vector weight m' = —1, 0, 1, as well as the body label n = 1, 2 of the spins; n = refers to quantities entering 
the spin-free part of h for which we also set m' = 0. As a result, for precessing binaries, the integral to compute him 
takes the form: 

1 2 „ 
m' = — 1 n=l 

where Xq^ = 1, Xo,„' = Xn' (^r n' = 1,2), A:_i^„/ = Q,, Xi,„/ = Cn', Ao,o ^o, Kqj^' /i^,, X-i,™' = fc„' 
and i^i^n' = fc*/. By means of the change of variable (j) (j) + a, we are able to factor out a complex exponential 
g-i(in+m )a contain all the dependence in a. Let us now focus henceforth on the case where the waveform has 

been expanded in powers of l. As we shall explicitly see below, [see Eq. (|4.16p ]. the /i£m's are then made of: (i) a 
spin-free piece proportional to e"""", (ii) two spin pieces proportional to e~'(™~i'" and to (* or Q respectively, (iii) 
two spin pieces proportional to e~'(™+i)" and to Ci or (2 respectively, (iv) two spin pieces proportional to e~"^" and 
to Xi or xi respectively. In contrast to what happens in the non-spinning case, him is not in general proportional to 
g-im* gxcept for the terms that are free of since both n and A reduce to trigonometric functions of $ -I- a as t ^ 0. 
The contributions to the polarization modes that are linear in l involve couplings of the type )(*+") cos$ or 

g-i(m+m )(*+q) sjj^<j) because the terms of first order in l entering n and A can only be linear combinations of icos$ 
or tsin(f> (or cquivalently ie***). Couplings like e^''-™^™ cos" $ sin'' with a, 6 G N, arise at higher orders 
making the dependence in <I> more complicated. A close inspection to the results below [see Eqs. (|4.16ap - ()4.16r|) with 
^' — > $] confirm these expectations. Beware that our mode normalization is tuned to factor out the exponential 
factors e~*™". 

A useful property of h coming from the arbitrariness of the body labeling is that it must be invariant in the 
exchange of particles 1 and 2: mi 7712, xi X2, n ^ — n, v —>■ —v. Under this transformation, the direction of 
the angular momentum L remains invariant hence A —A. The orbital frequency o^orb = (v • A) is unchanged as 
well as the direction of the total angular momentum, due to its structure and parity. Therefore, the phase $ becomes 
(xl, — n) = <f> + TT whereas the angles a and l are unaffected. This yields the relation 

h{0, cf,, c, a, $, xf,xlxl Xl Xlx^) = [h{e, 0, a, $ + tt, xl. xl.xl xl x'2,xl)]s^-s ■ (4.12) 

The previous identity may be supplemented by another one which originates from the classical parity invariance of 
physics: for any given time instant t, the waveform resulting from the stress-energy tensor parametrized by the world- 
lines x„ and the spins Xn must have the same value at point x as the waveform resulting from — x„ and +Xn at point 
— x. Taking into account the transformation of the polarization vectors under parity, this means for the function h: 

h{e, 0, t, a, xl,xl, Xn) = r(7r - 61, + TT, t, a, $ + tt, xl,xl, Xn) ■ (4-13) 

The above formula allows us to express the modes hgm in terms of the modes hi-m by performing the change of 
variable 6 n — d and (f) (p+ir in Eq. (|4.6p . The first factor of the integrand can be then rewritten as h {9, 0, t, a, $-f 
7r,x^,Xra)Xn) making use of Eq. (|4.13p . The second factor _2F^'"*(7r — 6*, </)-!- tt) may be transformed by means of two 
important symmetry properties of the spin- weighted spherical harmonics: _2y^'"(7r — O^ip + tt) = (— l)^+2y^™(^', '/') 
and +2Y^"'{e,(j)) = (-1)™_2F^-'"*(6I,0), which leads to the new expression (-l)™+^_2r^ "'" (6', (/)) for this factor. 
As a consequence, the link between him{^) = J dD, h{6, 0, l, a, $, Xn, Xn: Xn) -2Y^"^*{d, (f>) and /i^ _„($) is given by 

hirrm = (-1)"+^ J dnh*ie, Cjy, i, a, $ + TT, Xn, Xl, Xl) -2^' -"(0, 0) = (- l)^+"/l| ($ ^ ^) (4.14) 



The explicit expressions for the modes hem are obtained by inserting Eq. (|4.9p into the surface integral 
normalize them in such a way that the leading order mode starts with coefficient 1. Posing 



Dl \ 5 

and expanding in the inclination angle t, we arrive at ^ 
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{x: " *X^ + <^(x: - *xD) + e-'("+*) (x^ + ^Xl + SiXs + ^Xl)) 



— vt(5sin^' + 2i^cos2* - — v^tsin^'(xa + 5x1) ( + ^\ ^] > 
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12V14 I 



SVL 



+ 

11 135 



^<5(4 + i.)-16e-("+*)Kx:+*X^) 
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(5ut cos * + O 



2\/42 



^/f^(l-3.)+o(^|, 



8/10 



|^.3,(_i + 2.)+.Ml-3.)(e--^ 



^.2(l-3.)+0(- 



63 
21 V 2 



6v^ 
20 



(-l + 2i/)+w2i(l-3i.)e" 



(4.16a) 

*xD 

(4.16b) 

(4.16c) 

(4.16d) 
(4.16e) 

(4.16f) 
(4.16g) 
(4.16h) 

(4.16i) 

(4.16j) 
(4.16k) 

(4.161) 



In the case of spins aligned or anti-aligncd with the Newtonian angular momentum, the modes (2, 2), (2, 1) and (3, 2) were also computed 
in Ref. |70|| . We fully agree with their results. 
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.55 = -^^^5v^{l -2u)+o(\] , (4.16m) 



96V66 

hi = O , (4.16n) 



h2 = ^ (^) ' (4-16p) 
hi = -—^^j=5v''{\ ~2v)+0 , (4.16q) 



288V385 

ho ^ (^) • (4-16r) 

In Appendix [B] wc display the modes /122, /133, and /i2i for generic inchnation angle t. The modes hira for m < 
are derived from Eq. (|4.16p by means of the relation /i^ _„i($) = (— l)^/i|^($ + tt). The non-precessing expressions 
are obtained by setting u = Q and a = tt. Notice that when him{^) does not depend on we have simply 
hi-rn = (~l)^^|m- "'^^'^ Comparison with the modes of Refs. [13, HI] in the non-spinning case, it is important to be 
aware that the origin of the azimuthal angle there differs from ours by — 7r/2, which produces an extra factor (— «)™ 
(respectively i™) with respect to us in the modes (respectively in the spin- weighted spherical harmonics). 

Finally, let us emphasize that the (£, m) modes defined by Eq. (|4.5p depend on the particular choice of the source 
frame. In fact, they are functions of the spin and angular momentum components with respect to the (ic, y, z) basis 
introduced in Sec. [Ill As there is no canonical way to fix the reference frame for precessing binaries because of the 
secular but perpetual variation of the direction J/|J|, it is important to be able to relate the hi^s given in Eq. (|4.16p 
to the polarization modes h'^^ computed in another frame with different polarization vectors. Under a passive rotation 

(cosA -sinA 0\ / cosB sinS \ / cosF -sinF 
sinA cosA 10 sinF cosF 

1/ y-sinS cosB / V 1 

the (i,rn) modes transform in the same way as they would in the case of a standard spherical harmonics decomposition 
[tsI [79}. In fact, the spin- weighted —2 spherical harmonics are precisely devised to ensure this property for the modes 
of a spin-weighted —2 object [tI]. The law of transformation for the /i^m's is given by 

e 

h'i,„{^\a\i',x':',x/,x:')^ E An™'(AS,r) /.,„.(<&, a, i,x?^,xLxD, (4-17) 

rn=-e 

where the primed quantities refer to the new frame and where -D^/,„ is the unitary Wigncr matrix (76j 



D':^,^{A,B,r) = {-iy-'^-^_„,Y'"\B,A)e""'^ (4.18) 
with the convention of Landau-Lifchitz [tJ ■ The new angles read 

= arccos cos _B cos /, — cos(F + a) sini? sin t , (4.19a) 

cos A r cos L sin B + cos B cos(F + a) sin il — sin A sin t sin(F + a) 
a = arccos cos a = arccos ^ 

yl — ( cos B cos L — cos(F + a) sin B sin i) 

if cos t sin A sini? + sin t [ cos B cos(F + a) sin A + cos Asin(r + a)] >0, (4.19b) 

a' = 27r — arccos cos a' otherwise, (4.19c) 
$' = arccos cos $' 

cos B cos $ sin t — sin B sin F (cos t cos $ sin a + cos a sin <!>) + cos F sin B (cos t cos a cos $ — sin a sin $) 

1 — ( cos B cos i — cos(F + a) sin B sin i) ^ 

if cos$ sini?sin(F + a) + (cos t cos(F -|- a) sin B + cos i? sin i) sin $ > , (4.19d) 
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FIG. 2: The left panel shows the inclination angle of the orbital angular momentum relative to the total angular momentum, 
L, as a function of 2A/"i> for binaries with mass ratios 1:1 and 4:1, having initial spin orientations relative to the orbital angular 
momentum: SpiuA = {^i = 7r/2,(/)i =0,02 ~ 7r/2, (j!)2 = 7r/2} and Spine = {^i = n/6,(t>i = 7r/4, 6'2 = 7r/6, 02 = tt}- The right 
panel compares the modulus of h22 for the two precessing spin configurations with the non-spinning, aligned and anti-aligned 
cases for equal masses. The computations use waveforms accurate to 1.5PN in amplitude and phase evolved with the precession 
equations at 1.5PN order [see Eqs. ([3TT|| . (f3T8|l . (|T20|l . and (fMI]) ]. Note that these plots (and those of Figs. 3 and 4) begin 
at 2 Af $ = 0.02 which is approximately where the dominant second harmonic from a binary of total mass 16 Mq enters the 
LIGO band at 40 Hz and where the second harmonic from a binary of total mass 6.5 x 10^ Mq enters the LISA band at 10^* 
Hz. 



$' = 27r — arccoscos$' otherwise. (4.19e) 

When the direction of the total angular momentum used to built the new frame coincides with that of Jq, which 
results in the equality z' = z, the Euler angle B vanishes. Then, it can be checked from Eqs. (|4.19dp and (|4.19ep that 
<&' = <I> as expected. 



V. FEATURES OF GRAVITATIONAL- WAVE MODES FOR PRECESSING BINARIES ON NEARLY 

CIRCULAR ORBITS 

We now study how spin effects change the waveform modes for generic precessing binaries. We consider two 
maximally spinning configurations with mass ratios 1:1 and 4:1. We label the spin configurations with the angles 
{6*1, 411,02, 412}, where {6*;, (f>i} describe the orientation of the spin vector of the i^^ body relative to the orbital angular 
momentum in the initial configuration, which we take to be a circular orbit with Mworb = 0.001. We use the 
full expressions for the /ifm's (i.e., the expressions that have not been expanded in i) as given in Appendix IB] and 
normalized following Eq. (|4.15p . but we replace 5* with 

After evolving through 1.5PN order all dynamical quantities they depend on [sec Eqs. (|3.1ip . (|3.18p . p.20p . and 
(|3.21|) ]. we compute the modulus — more often referred to as the absolute value — of a sample of modes. Considering 
the complicated structure of the hg^'s, their qualitative behavior in the presence of spins is discussed here in terms of 
the t-expanded formulae (j4.16p . Let us focus on two spin configurations. The configuration SpiuA = {7r/2, 0, 7r/2, tt/2} 
has both spin vectors in the orbital plane, meaning a relatively large inclination angle. The configuration Spine = 
{tt/B, 7r/4, 7r/6, tt} has a smaller component of total spin transverse to the orbital angular momentum, hence a smaller 
inclination angle. In Figs. [3] and [3] we plot the amplitude of the hirn over the frequency range 2M<i> = 0.02-0.15, 
the upper frequency being reached roughly 2 cycles before merger, for an equal-mass, non-spinning binary [63 |. For 
a binary of total mass WMq, the dominant second harmonic varies over the frequency range 40-300 Hz. For a 
6.5 X 10^ Mq binary, this range is shifted to 10"^-7.5 x lO"'* Hz. 

The left panel of Fig. [5] shows the inclination angle t as a function of the dimcnsionlcss frequency 2il/(f>. We see 
that the inclinations are much larger in the case of 4:1 mass-ratio than for equal masses. Moreover, the inclination 
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FIG. 3: We plot the modulus of the ^ = 2 modes for mass ratios 1:1 (left panel) and 4:1 (right panel) with the spin configurations 
described in Fig. [H The computations use waveforms accurate to 1.5PN order in amplitude and phase evolved with precession 
equations at 1.5PN order. The dashed lines are the larger t configuration (Spin^) and the solid lines are the smaller t 
configuration (SpiUg). We see that as t increases, the modulus of /i22 decreases while the modulus of the other 1 = 2 modes 
increases. This effect becomes more pronounced when the mass ratio is more extreme. 



increases monotonieally in the equal-mass case, whereas the 4:1 mass ratio exhibits nutation, since the inclination 
oscillates, but grows on average. These observations can be explained as follows. At early times, when the binary 
has a large orbital separation, we have |L| ^ |S|, where S = Si + S2 is the total spin, so that J = L + S and L are 
nearly aligned. Radiation reaction causes |L| to decrease, making J move away from L and toward S. This is why the 
inclination angle, t, grows on average as the frequency increases. The absence of oscillations for the inclination angle 
i in the equal-mass case can be explained by the fact that we are evolving the dynamics, in particular the precession 
equations (|3.18|) . through 1.5PN order, i.e., wc arc neglecting spin-spin effects. Due to the equality fJi = ^2 at this 
accuracy level, the precession equations simplify then to a single equation of the form dS/dt — fl x S. In the absence 
of radiation reaction, S precesses around a fixed direction with a constant frequency, and the inclination is constant 
(apart from the increase produced by radiation reaction). For unequal masses this symmetry does not exist, with the 
consequence that one must solve two coupled equations for Si, S2 instead of a single equation for S. The motion of 
the spin vectors is thus more complicated. Schematically, they rotate about a fixed direction while also bobbing up 
and down (soj . 

The right panel of Fig. [2] plots, for the case of equal masses, the absolute value of the h22 mode normalized to 
its Newtonian order expression, h22, for both precessing spin configurations as well as the non-spinning, aligned and 
anti-aligned cases. One interesting feature is that the aligned and anti-aligned cases do not bound the absolute 
value of the modes for generic spin configurations. This and other features of the plot can be understood from 
the L expansion (|4.16ap of /i22, which contains four spin corrections. The first correction, —t? 12. is zero for the 
aligned and anti-aligned cases, while it decreases the absolute value of the h22 mode for all other spin configurations. 
If L is comparable to 1 radian, it can be a significant correction. The second correction, (l/3)e**i5?;t, vanishes 
for equal masses. For unequal masses, it interferes with the non-spinning terms and, because it has a different 
dependence on the orbital phase, produces oscillations in the absolute value of /i22- Next, the IPN order spin 
correction, — (i;^/2)e*(""'"*^[xa — + i5(x? ~ *Xs)]i generates oscillations that depend on the spin vector components 
transverse to the total angular momentum. Finally, the 1.5PN order spin correction, 45xa/3 + (4/3)(— 1 + v)xi\^ 
lowers (raises) the absolute value of /i22 for spins aligned (anti-aligned) with the total angular momentum. It is solely 
responsible for the spread between the aligned and anti-aligned cases, as the other corrections all vanish then. 

A similar analysis can be applied to understand the behavior of the other modes. As an illustration, we plot in Fig. [3] 
all of the 1 = 2 modes for mass ratios 1:1 and 4:1. The /121 mode ()4.16bp is zero for non-spinning equal mass binaries. 
However, it contains several spin corrections and can have significant amplitude for precessing binaries, particularly 
for large l. It can exhibit complicated modulation, as its different spin corrections interfere with one another. The 
/i2o mode (|4.16cP also has several spin corrections, most notably 2 i? cos\I'. This correction is responsible for the large 
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FIG. 4: We plot the modulus of the 1 = 3 modes (left panel) and £ = 4 modes (right panel) for equal masses with the spin 
configurations described in Fig. [5] The computations use waveforms accurate to 1.5PN order in both amplitude and phase 
evolved by means of the 1.5PN precession equations. The dashed lines refer to the larger t configuration (Spin^), the solid lines 
to the smaller iota configuration (Spiug). We see a redistribution of power among the modes similar to Fig. [3] As t increases, 
the largest modes for the non-precessing cases (1^32 1 and |/i44|) become smaller, while the other modes become larger. 



oscillations in the absolute value of /i20- Note that in the late stages of the inspiral evolution for the 4:1 mass-ratio 
case, where t ~ 1 radian, these oscillations in the absolute value of peak near the absolute value of /122. The other 
spin corrections in /120 are responsible for the further modulations of the absolute value. 

Figure 0] plots the absolute value of the £ = 3 and £ = 4 modes for equal masses. Note that these modes are about 
two orders of magnitude smaller than the £ = 2 modes. This remains true in the non-spinning case for the mass ratios 
we consider. In the case of non-spinning, equal mass binaries, for £ = 3 only the /132 mode (|4.16eP is non-zero. The 
1.5PN order spin correction decreases (increases) this mode's absolute value if the spins are aligned (anti-aligned) 
with the total angular momentum, in a similar way as in the mode. For unequal masses, the /132 mode (|4.16eP also 
has an interference term proportional to l. The other £ = 3 modes are non-zero for generic precessing binaries, and 
generally have larger absolute values for larger inclinations. In the left panel of Fig. IH we indeed observe that |/i33|, 
|ft,3i|, and |/i3o| are greater for the configuration Spin^ than for Spiug. For £ ~ 4, only the /i44 and ^,42 modes (j4.16hp . 
( |4.16j[ ) arc non-zero for non-spinning, equal-mass binaries, with /144 being the largest. Though the t-expanded form 
of the /i44, /142 and ft.40 modes in Eq. (|4.16p do not have any spin corrections through ©(l/c"^), we do see spin effects 
in all of the £ = 4 modes when we plot the full expressions accurate through . The v'^ and coefficients in the full 
expressions for the /144, /142 and /140 modes depend on the inclination l but not on the spin vector components and 
this dependence is such that if we treat /, as a 1/c correction by performing a Taylor expansion in powers of l, then the 
spin terms are proportional to v'^ i^, l and higher order in l. They are thus considered as higher order corrections 
in the l expansion, though they are present when we expand only in powers of v. None of the £ = 4 modes contain 
any spin corrections proportional to the spin vector components through order v^. Those corrections would appear 
only at higher order in v. 

While wc have plotted the £ = 3 and £ — 4 modes solely for equal masses, we have also studied these modes for 
4:1 mass-ratio binaries. We find that they are affected by the change of mass ratio in much the same way as the 
£ = 2 modes: the redistribution of signal among the £ = 3 and £ = 4 modes is more pronounced for asymmetric 
binaries than for equal mass binaries. However, even for 4:1 binaries with large l (SpiuA) spin configuration, all of 
the i = 3 and £ = 4 modes are still one or two orders of magnitude smaller than /122, whereas the /i2i and /i2o modes 
can be comparable to /i22- The reason is essentially that /i2i and /i2o have l corrections at leading order in v, while l 
corrections to the £ = 3 and i = 4 modes all appear at higher order, and so they do not have as strong an effect as 
for /121 and /i2o- 

By examining the absolute values of the h/m modes for precessing binaries, we see that they can be significantly 
altered by the motion of the orbital plane relative to the frame used to perform the mode decomposition and the 
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signal may be redistributed among the modes as observed in Ref. [79|. This suggests that all modes, not just the 
dominant ones for non-spinning binaries, are needed to accurately describe the waveforms emitted from precessing 
binary systems, especially for asymmetric binaries and binaries with large inclinations where this redistribution of 
signal among the modes is most dramatic. 

Let us close this section with a few comments about the applicability of the full and t-expanded expressions for the 
modes and the possibility of combining them with higher order non-spinning corrections. We find that the absolute 
values of the full and t-expanded him modes are often quite close to each other for relatively small t. For inclinations 
less than half a radian (30°), the difference in I/122I is typically of a few percent. For t comparable to or larger than 
a radian (60°), significant differences between the full and t-expanded modes develop and the absolute values may 
differ by ~ 10-100% when l > 1 radian. Nonetheless, the t-cxpandcd modes are very useful in understanding the 
qualitative behavior of precessing binaries, even for inclinations ^ 1 radian, albeit they should not be used for precise 
quantitative studies of binaries with large inclination angles. 

In Refs. [27I . [H], expressions of the modes arc given to 3PN order for non-spinning binaries. Wc have compared 
their absolute values to that of the corresponding quantities truncated at 1.5PN order, which has shown us that they 
typically differ by ~ 1-10%. For example, the absolute value of the 1.5PN and 3PN order modes for equal-mass 
binaries differ by less than 1% at 2M $ = 0.01 and by about 3% at 2M$ = 0.05 or 2M$ = 0.12. The other modes 
typically have a larger difference. The 1.5PN order and highest known order absolute values for /132, /144 and /142 differ 
by about 5-15% over this same frequency range for equal mass binaries. For 4:1 mass ratio binaries, the differences 
in absolute value arc similar. 

Known higher-order non-spinning terms can actually be included in the amplitude if enough care is taken. In 
the non-prcccssing case, the modes h^m are proportional to e^""*, as in the non-spinning case, because ii and A 
appearing in the strain tensor hij are trigonometric functions of the orbital phase $. In constructing the /i^m's from 
hij, they generate an exponential dependence on multiples of <&. However, for the case of a precessing binary with 
small t, n and A are trigonometric functions of $ -I- a. Thus, the /i^m's contain then all of the non-spinning terms, 
but with the substitution 'if ^ + a. The situation is different for the general precessing case. The vectors n and 
A depend on <I>, a, and l, and the resulting have a complicated dependence on all three of these quantities that 
cannot be simply related to the non-spinning case. These considerations show that it is only for binaries with a 
small inclination (or no inclination) that we can readily construct the himS with spin effects up to 1.5PN order and 
non-spinning correction s up to 3PN order. For spins (anti-) aligned it is trivial to add the higher-order non-spinning 
corrections of Refs. [l^, [Sj to the hirn given in Eq. ()4.16p . For precessing binaries with small inclinations, they can 
be added to our expressions in Eq. (|4.16p with the substitution 5* + a. For general precessing binaries, it is not 
so simple to include higher-order non-spinning corrections to the full expressions for the him given in Appendix |B] To 
do this properly, we would need the spin terms at the same order as the non-spinning terms and repeat the derivation 
of the hi„i to a higher order. 

VI. READY-TO-USE FREQUENCY-DOMAIN TEMPLATES FOR SPINNING, NON-PRECESSING 

BINARIES 

A. Gravitational-wave polarizations in time domain 

In the non-prcccssing case, the orbital angular momentum points in a fixed direction which we take to be the z-axis 
(see Fig. [T]) and the spins are either aligned or anti- aligned with it. The basis vectors of the orbital plane, Xi,(t) and 
yi(t), are constant in time. They can be freely chosen to be any pair of orthogonal unit vectors in the x-y plane. 
Here, following the convention of Ref. [2^, we choose = P = N x Ln/|N x Ln|, so that the phase is zero at the 
ascending node (where the orbital separation vector crosses the plane of the sky from below). This is equivalent to 
setting t = and a = tt in Eqs. (|3.12p . (jA2p . and (|A3[) . Note also that since the orbital plane remains fixed, the phase 
$ defined through Eq. (|3.10p coincides with the standard definition of the orbital phase, that is 

Worb-$. (6.1) 

In the non-precessing case, the vectors in terms of which the GW polarizations are expressed originally take the 
simpler form 

fi = (sin$, -cos$, 0) , (6.2) 

A= (cos$, sin$, 0) , (6.3) 

N = (sin 6*, 0, cose*) , (6.4) 
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(0,0, 1) 



(6.5) 



By plugging the expressions (|6.2p - (|6.5p into Eq. (|3.3[) and taking the combinations given in Eq. (|2.ip . we obtain an 
equation similar to Eq. ()3.12p . The spin-dependent IPN, 1.5PN and 2PN order polarization coefhcients read 
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(6.6a) 
(6.6b) 

(6.6c) 
(6.6d) 

(6.6e) 
(6.6f) 



where SS labels the spin(l)-spin(2) contributions. ^ In the equations above, we use the shorthand cg = cos9 and 
Sg = sin0. Note also that the phase \1/ is the shifted orbital phase that relates to $ at our accuracy level as 



= $ 
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(6.6g) 



where can be chosen arbitrarily. Expressed in terms of the orbital phase <I>, the GW polarizations would contain 
terms logarithmic in tJorb, arising from the propagation of the tails. However, introducing the phase ( |6.6gP , they are 
all absorbed (up to the 2.5PN order we arc considering) into the phase variable [25| . 

The spin-dependent polarizations (|6.6ap - (|6.6f|) were derived in Rcf. [s^ [see Eqs. (F24a)-(F25c) in that paper], 
although the 1.5PN and 2PN order cross polarizations had an erroneous sign, which is corrected here. 

B. Spin-orbit effects at 1.5PN order and spin-spin effects at 2PN order in the frequency-domain 

gravitational-wave amplitude 



Writing h{t) = + and collecting terms by PN order and by sines or cosines of harmonics of the orbital 

frequency, we can write the time-domain strain in the compact form: 
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J2^k^" K"^ cos(fc*(t))+/3f^ sin(A:^(0)] + 0{v^) , 

n=0 k=l 

E E ^k^" cos(fc ^{t)) + cos(fc ^{t) - ^)1 + Oiv') , 



n=0 k=l 



(6.7) 



where n/2 is the PN order and k labels the harmonics of the orbital phase. The PN expansion parameter is defined 
as = {2T:MFy/^, with F = Worb/(27r). We shall denote the GW frequency by /. For the A;'^ harmonic, we have 
then the relation f = fk ^ kF, so that 



2 TT A/ ^ 



(6.8) 



Given a function of the form h{t) = A{t) cos (f>{t), where (j){t) is a monotonically increasing hmction satisfying 
d\iiA{t) /dt <^ d(j}{t)/dt^ we can compute its Fourier transform by applying the stationary-phase approximation (SPA): 



1 



27r 



h{f) ~ -A{t{f)) 



,i(2 7r/t(/)-0(t(/))-7r/4) 



(6.9) 



We remind that spin(l)-spin(l) and spin(2)-spin(2) effects in the waveform polarizations are currently miknown. 
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t{f) being defined here for each frequency / as the value of t for which {d(j)/dt)(t) = 27r/. In a similar manner, we 
apply the SPA to each term in the sum of Eq. (|6.7p . Moreover, for each harmonic of the orbital phase, we expand 
the factor inversely proportional to the second time derivative of the orbital phase entering Eq. (|6.9p in a PN series 
of the form 




k^) = \l^ [1 + ^2 Vi + 53 Vi + 54 + S, + 0{V^)\ , (6.10) 



with 
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5. = .f-iI^ + ^.V (6.11) 

The expansion (|6.10p without spin corrections in the amplitude was first given in Ref. [l^l . We have added to it the 
leading order SO corrections through 1.5PN order and the spin(l)-spin(2) SS corrections appearing at 2PN order. In 
principle, SO corrections at 2.5PN order and spin(l)-spin(l), spin(2)-spin(2) SS corrections at 2PN arising from the 
spin contribution to the orbital frequency are also present. However, when calculating spin terms in the frequency- 
domain amplitude, we neglect them because they have not been calculated yet beyond the 1.5PN order in the time 
domain amplitude. The spin contribution at 2PN and 2.5PN order to the Fourier domain amplitude is not complete 
unless we take both into account. 

Defining the frequency-dependent SPA phase as 

*SPA(./) = 27r/t(/) -*(/), (6.12) 

the frequency domain waveform with amplitude corrections containing SO effects through 1.5PN order and spin(l)- 
spin(2) effects through 2PN order is 

«=0 fc=l ^ ^ 



A/2 I5nu 
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where the index n denotes the PN order and the index k the harmonics. Explicit expressions for the C^" can be found 
in Appendix |D] The non-spinning terms in the amplitude agree with Ref. [46j , although we have written them in a 
different, more explicit manner. Notice that recently the non-spinning amplitude corrections were calculated through 
3PN order [23|, but in this paper we restricted the computation to 2.5PN order. 



C. Spin-orbit effects at 2.5PN order in the frequency-domain gravitational- wave phase 

For matched filtering, it is best to know the GW phasing at the highest PN order. We now derive the SO 
contributions to the SPA phase through 2.5PN and the SS contributions {including spin(l)-spin(l) and spin(2)-spin(2) 
contributions) to the SPA phase through 2PN order. 
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The PN expansion of the SPA phase ^spa{F) can be obtained from the PN expansions of the binary center-of-mass 
energy, E, and GW flux, F, via the energy balance equation 

Using d^/dt = 2 7rF = v^/M, we can re-write the energy balance equation as the differential equations 

dt^-^^dv, (6.16) 
dv J- 

and 

d* = _^l!^d.. (6.17) 

dv T M ^ ' 

The quantities E and T are known as power series in v = {2tt M F)^/^ . The non-spinning terms in the expansions 



of E and T have been calculated by Refs. |2£ . |81| . |82|. l83l. l84l. 18511 . while the spin contributions to these quantities 
through 2.5PN order were derived by Refs. 22 



l23l l30l . |3J, |35|, |7J] . The center-of-mass energy and the flux read 



where the coefflcients Ei and Fi are explicitly given in Appendix [Cl By inserting Eqs. (|6.18p . (|6.19p into Eqs. (|6.16p . 
(|6.17p . we obtain rational function approximations to the integrands. We then find the Taylor series of the rational 
functions and integrate up to some reference frequency, often chosen to be the time of coalescence, when the orbital 
frequency formally diverges. Thus, we obtain PN approximations of the form 

t{v) = tc- (2v-' + j2t,v=A , (6.20) 
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(6.21) 

where the tj coefficients are linear combinations of products of the Ei and Fi. Plugging Eqs. (|6.20p . (|6.21[) into 
Eq. (|6.12p . we obtain the following expression for the SPA phase through 2.5PN order ^ 

^spa{F) - 2^Fi, - + — (2^XF)-5/3 J 1 ^ / __ + \ {2nMFf/^ + (4/?- 16^) {2^MF) 
15293365 27145 3085 , \ , / 38645 6£ 



508032 504 72 



\ , ,4/., /38645 65 \ , , , 

-lOaj {2T:MFf/^ +{-^'k-—t:v \ (1 + 3 log(z;)) x 



(27rMF)5/3| ^ (g 22) 

where the 1.5PN SO phase corrections are contained in [3, the 2PN SS corrections are contained in cr, and the 2.5PN SO 
corrections are contained in 7. Note that (3 and the spin(l)-spin(2) contributions to a were previously known [l3l.[l^. 
while we have calculated the spin(l)-spin(l) and spin(2)-spin(2) contributions to a and the 2.5PN SO corrections to 
the SPA phase using the results for the center-of-mass energy and GW flux of Refs. [1^, [13, [13, [zj • Explicitly, these 
corrections are 

/3 = f^-yH X.-LN + ^<5Xa-LN, (6.23) 



The non-spinning terms in the SPA phase through 3.5PN order can be found in Ref. [86 
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We note that these expressions are only valid when both component spins are ahgned or anti-ahgned with the orbital 
angular momentum. The spin(l)-spin(l) and spin(2)-spin(2) contributions to a were also derived in Refs. [3^.[33j and 
we found full agreement with them. 



D. Features of frequency-domain non-precessing waveforms with higher hcU-monics 

We now discuss some interesting features of the spinning, non-precessing waveforms derived in Sec. IVIBI Several 
papers have studied the effect of higher harmonics in the amplitude corrections of non-spinning binaries observable 
by ground- and space-based detectors [iel. IItI. Isills^. IssLIstI . 

One important feature of the higher harmonics in the waveform amplitude is that they can increase the mass reach 
of a detector [5l|. This is because high-mass binaries whose dominant second harmonic is not in the detector's sensitive 
band can have higher harmonics in band and therefore become visible to the detector. A closer look at Eq. ()6.14|) 
and Appendix ID] shows that spin corrections through 2PN order appear only in the first and second harmonics. In 
particular, the only SPA amplitude coefficients with spin dependence are Cf\ Cf^ and 4"^^ given in Appendix |D] 
Thus, in the non-precessing case spin corrections through 2PN order in the waveform amplitude do not affect the 
mass reach of the detector, and only affect binaries whose second harmonic appears in band. 

Another general feature of the higher harmonics is that they interfere with one another, typically destructively (46l . 
[Sll, HI]. For binaries that would be visible with Newtonian waveforms, this effect tends to decrease the signal to noise 
ratio (SNR). As we shall study in detail in this section, spin effects can play a role in this interference, either raising 
or lowering the SNR depending on the spin orientations. 

We define the power spectrum, P{f), as 



and the optimal SNR, p, as 



where fs is the low frequency seismic cutoff of the detector, and the upper frequency cutoff is taken to be the highest 
harmonic of the orbital frequency at the last stable orbit (LSO) which for simplicity wc choose to be the LSO of a 
test particle in Schwarzschild, 

Note that the A:*'^ harmonic ends at /c^lso as enforced by a step function 9{kFi^so ~ /) [see Eqs. (|Dip - (|D21ll in 
Appendix ID] . In Eqs. ()6.26p . (|6.27p . we denote with Sn{f) the noise power spectral density of the detector. For 
Advanced LIGO, we take the spectral density to be Eq. (4.3) of Ref. [43| and fix fs = 20 Hz. For LISA, we use the 
so-called effective non-sky-avcragcd spectral density given in Eqs. (2.28)-(2.32) of Ref. [s^. We do not consider the 
orbital motion of the LISA spacecraft [s^ and consider only the single detector configuration In the presence of 



It should be noted that in our model, though wc do not perform an average over the antenna pattern functions, wo do not account 
for the orbital motion of LISA either. In this sense, our model falls in between the pattern averaged and non-pattern averaged cases 
described in Ref. ISSlI 
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FIG. 5: We compare the power spectra computed with the Newtonian amplitude waveform (red dashed line) and the 2.5PN 
waveform with 1.5PN SO and 2PN SS effects included (blue, continuous line). In the left panel we consider a typical source 
for LISA, a binary with total mass (10® + 10'')Mq, and spins maximal and aligned with the orbital angular momentum. In the 
right panel we consider a typical source for Advanced LIGO, a binary of total mass (30 + 30)Mq with spins xi = 1, X2 = 0.5 
aligned with the orbital angular momentum. Note that the k*^ harmonic ends at kFhso, and these frequencies are marked 
by the vertical dashed lines on the graph. The spectrum of the 2.5PN waveform is much simpler in the equal-mass case than 
unequal mass case because in the former case all non-spinning odd harmonics are suppressed. 



higher harmonics, the lower and upper cut-off frequencies are chosen following Sec. lllA of Rcf. For LISA we 

assume an observation time of one year, and the orbital frequency at the beginning of observation to be Eq. (3.3) of 
Ref. [5l|. As explained in Ref. [5l|, this can be implemented by multiplying the k^^ harmonic by the step function 
d{f — kFin) where Fm is the orbital frequency at the beginning of observation. Finally, because of the 60° angle 
between EISA's arms, we use h{f) (\/3/2) h{f) in Eqs. ([07|) in the case of LISA. 

All tables and figures in this section, refer to a binary with orbital angular momentum inclined relative to the line 
of sight by = 7r/3, sky location 9 ^ <f> ^ tt/6 and polarization angle -0 = 7r/4 [see Eqs. (|2.4p . (|2.5p ]. We have verified, 
by considering random values for the four angles, that the qualitative trends reported in this section are generic and 
do not depend on the specific values of them (see a detailed discussion at the end of this section). Regardless of the 
PN order of the amplitude, all waveforms use the SPA phase with non-spinning terms up to 3.5PN order [s^, and 
spin terms up to 2.5PN order, as given in Eqs. (|6.22p - (|6.25p . In the case of Advanced LIGO (LISA) we consider 
binaries at a distance of 100 Mpc (3 Gpc). Moreover, all masses and distances refer to the redshifted quantities. 

The destructive interference between different harmonics can be seen in Fig. [51 The Newtonian waveform's power 
spectrum is simply proportional to f~''^^ / Sn{f)- The higher harmonics present in the 2.5PN waveform introduce 
oscillatory cross terms that on average lower the power. Notice that although the higher harmonics extend the 
observable frequency band significantly, the power beyond the cutoff of the second harmonic, being at a higher PN 
order, is suppressed by one or several orders of magnitude. These features explain why the SNR listed in Tables [H 
im tends to decrease as the PN order increases for the range of masses we consider. For equal-mass binaries, all 
non-spinning odd harmonic corrections arc suppressed because the latter are proportional to 5 which is zero for equal 
masses [see Eqs. (jDip - (|D2ip in Appendix [P] . This is not true of spin-dependent amplitude corrections. For example, 
the first harmonic has a spin dependent amplitude correction at IPN order which docs not vanish for equal mass 
systems unless spins are equal and aligned with one another [see Eq. (|6.6ap ]. 

Tables U |ll] show the SNR for the case of maximal spins both aligned or anti-aligned with the orbital angular 
momentum. From the bottom three rows of Tables |T1 |TT1 we see that, depending on the spin orientation, the 2.5PN 
amplitude corrections with spins can have SNR ~ 10% higher or lower than the 2.5PN amplitude corrections without 
spins. We caution that this ~ 10% change in the SNR from spin corrections is only meant as a bound on spin effects 
for spinning, non-precessing binaries. As we have seen in Sec. |Vl the affect of spin corrections on prcccssing binaries 
is not bounded by the cases of maximal spins aligned and anti-aligned with the orbital angular momentum. 

Quite interestingly, the 1.5PN SO and 2PN SS corrections arc far more important than the IPN SO correction in 



23 




FIG. 6: For a binary of total mass (30 + 30)A'[q with spins xi = 1, Xz = 0.5 aligned with the orbital angular momentum (the 
same binary of the right panel of Fig. [5}, we show the power spectra up to 2Flso. We plot the power spectrum for the waveform 
through 2.5 PN order with no spin corrections (cyan solid line) and with SO corrections through 1.5PN (that is, IPN and 
1.5PN) and SS corrections at 2PN order (dark blue solid line). We also plot power spectra for the waveform with Newtonian 
amplitude (red dashed line), Newtonian amplitude plus the IPN SO correction (black dotted line), Newtonian amplitude plus 
SO effects through 1.5PN (green, dot dashed line), and Newtonian amplitude plus SO corrections through 1.5PN and the 2PN 
SS correction (magenta, double-dot-dashed line). The 1.5PN SO and 2PN SS effects raise and lower the power in the dominant 
harmonic while the IPN SO effect merely changes the modulation pattern up to its cutoff frequency of -Flso. Vertical dashed 
lines mark the frequencies Flso and 2_fLSO. 



terms of their effect on the power spectrum and the SNR. Notice that in Tables [H [IT] the IPN SO term always has 
httle or no effect, while the 1.5PN SO term changes the SNR by 10%, and the 2PN SS term changes the SNR for 
the equal-mass binary. The reason the IPN SO term is less important is that the 1.5PN SO and 2PN SS terms are 
corrections to the second harmonic, so they increase or decrease the power in the dominant term. On the other hand, 
the IPN SO term is a correction to the first harmonic. Thus, it is merely a perturbation to the dominant signal, and 
only in the lowest part of the spectrum where the first harmonic is observable. This is illustrated in Fig. [51 where we 
plot the power spectrum as a function of frequency (up to 2Flso) for different spin contributions for the (30 4-30)Mo 
binary system. We see that the 1.5PN SO and 2PN SS corrections add or subtract their power coherently with the 
dominant second harmonic. Their net effect is to shift the power spectrum of the full waveform upward without 
changing its shape. On the other hand, the IPN SO correction, which is proportional to the (sine or cosine of) half 
the dominant harmonic, simply changes the modulation pattern of the full waveform up to i^LSO (37 Hz). It should 
however be noted that the structures in the power spectra could be more complicated for asymmetric systems where 
the non-spinning terms proportional to cos 5* and sin ^! arc not suppressed. 

The 1.5PN SO term is typically the most important of the spin terms. This term is linearly proportional to the 
spins of the two bodies, as can be seen in Eq. (jD8|) . If the spins are aligned with the orbital angular momentum it 
increases the SNR. If the spin terms are anti-aligned with the orbital angular momentum it decreases the SNR. If one 
spin is aligned with Ln and the other anti-aligned, the body with the greater spin Si = mf Xij which is typically the 
larger body, dominates. Thus, the large body dictates whether the SO effect increases or decreases the SNR, unless 
the spin of the smaller body is much greater than the spin of the large body. This is illustrated in Table IIVI where 
the mass ratio mi : m2 = 10 : 1. The spin of the larger body is aligned with Ln and tends to increase the SNR while 
the spin of the smaller body is anti-aligned with Ln and tends to decrease the SNR. For a spin ratio xi '■ X2 — ^ '■ ^ 
there is a large increase in SNR due the larger BH. For a spin ratio 1:10, the larger BH still dominates and we get a 
small increase in SNR. For the spin ratios of 1:100 and 1:1000, the smaller BH is now able to overcome the larger BH 
and produce a net decrease in the SNR. 

The 2PN SS term decreases the power spectrum and SNR when the component spins are aligned with one another, 
and increases the power spectrum and SNR when they are anti-aligned with one another. The 2PN SS term has a 
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Advanced LIGO SNR 





(50 + 5)Mq 


(30 + 30) Mq 




(1,1) 


(-1,-1) 


(1,1) 


(-1,-1) 


Newt 


76.4 


76.4 


131.1 


131.1 


0.5PN 


84.3 


86.5 


131.1 


131.1 


IPN 


73.7 


77.4 


116.9 


115.2 


IPN + IPN SO 


73.7 


77.6 


116.9 


115.2 


1.5PN 


68.9 


72.7 


116.9 


115.2 


1.5PN + 1.5PN SO 


79.4 


62.6 


134.3 


98.2 


2PN + 1.5PN SO 


75.5 


59.1 


123.2 


88.0 


2PN + 1.5PN SO + 2PN SS 


75.2 


58.8 


121.2 


86.1 


2.5PN 


63.4 


66.7 


104.5 


102.9 


2.5PN + 1.5PN SO 


73.5 


57.4 


121.7 


86.2 


2.5PN + 1.5PN SO + 2PN SS 


73.5 


57.4 


121.7 


86.2 



TABLE I: For several binary configurations observable by Advanced LIGO we fist the SNR as the PN order of the amplitude 
corrections is varied. In each column we show the component spins (xi • Ln, Xz • Ln). We include all non-spinning, spin-orbit 
(SO) and spin-spin (SS) corrections up to the orders given in the first column. For example, 2.5PN + 1.5PN SO -I- 2PN SS 
means we include non-spinning amplitude corrections from Newtonian to 2.5PN order, IPN and 1.5PN SO corrections, and 
the 2PN SS correction. Regardless of the PN order of the amplitude, we always use the SPA phase with non-spinning terms 
up to 3.5PN order, and spin terms up to 2.5PN order, as given in Eqs. (|6.22|| - (|6.25|l . The binary is at a distance of 100 Mpc 
with orbital angular momentum inclined relative to the line of sight by 6 = 7r/3, sky location 9 — (j) — tt/G and polarization 
angle tp = 7r/4 [see Eqs. plijl . lp^ll ]. 







LISA SNR 






(2 X 10" + W*)Mq 


(10"* + 1O®)M0 




(1,1) 


(-1,-1) 


(1,1) 


(-1,-1) 


Newt 


382.6 


382.6 


2764.5 


2764.5 


0.5PN 


599.2 


598.7 


2764.5 


2764.5 


IPN 


620.3 


619.9 


2515.9 


2456.4 


IPN + IPN SO 


620.1 


619.9 


2515.9 


2456.4 


1.5PN 


512.2 


510.9 


2515.9 


2456.4 


1.5PN + 1.5PN SO 


552.0 


477.8 


2882.2 


2105.7 


2PN -f 1.5PN SO 


524.0 


452.9 


2612.7 


1861.5 


2PN + 1.5PN SO + 2PN SS 


524.0 


452.9 


2569.7 


1823.0 


2.5PN 


479.7 


478.2 


2286.8 


2230.6 


2.5PN + 1.5PN SO 


516.8 


448.9 


2646.1 


1891.3 


2.5PN + 1.5PN SO + 2PN SS 


516.7 


448.8 


2603.7 


1853.9 



TABLE II: For several binary configurations observable by LISA we list the SNR as the PN order of the amplitude corrections 
is varied. In each column we show the component spins (xi ■ Ln,X2 ■ Ln). We include all non-spinning, spin-orbit (SO) and 
spin-spin (SS) corrections up to the orders given in the first column. Regardless of the PN order of the amplitude, we always use 
the SPA phase with non-spinning terms up to 3.5PN order, and spin terms up to 2.5PN order, as given in Eqs. (|6.22p - (|6.25fl . 
The binary is at a distance of 3 Gpc with the same orientation as in Table U The binary masses and distances refer to the 
redshifted quantities. 



greater effect on the SNR and power spectrum than the IPN SO term, but is less important than the 1.5PN SO term. 
This is because it is suppressed relative to the 1.5PN SO term by a factor oi v/c and it is quadratic in the spins and 
proportional to the symmetric mass-ratio v. Thus, the 2PN SS term are most important for binaries with two large 
component spins and comparable masses. From Tables H Hill wc can see that the 2PN SS term has little or no effect 
on binaries with a mass ratio greater than 10:1. In Tables Hlllll for the columns with equal masses and spins aUgncd 
with one another, the 2PN SS term decreases the SNR by a few percent. For the binary in Table Hill we sec that 
the 2PN SS term increases the SNR by an amount comparable to the SO terms when the spins are maximal. As we 
decrease the spin magnitude, the SS effect is suppressed faster than the SO effect because it is quadratic in the spins 
while the SO effect is linear. 

Before ending this section we study how different values of the source position and incUnation angle can affect the 
SNR trends shown in Table H For (5 + 50)Mq and (30 + 30)Mo systems we calculated the SNRs at different PN 
orders in ampUtude for various random realizations of 0, </), "0 and 9 and for the spinning and non-spinning cases. For 
the spinning cases, when all the known spin effects are included at different PN orders, the trends across different 
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Advanced LIGO SNR 






(60 + 4O)M0 






(1,-1) 


(0.8, -0.8) 


(0.5, -0.5) 


(0.2,-0.2) 


2.5PN 


79.6 


81.9 


81.6 


82.0 


2.5PN + 1.5PN SO 


83.1 


84.6 


83.4 


82.7 


2.5PN + 1.5PN SO + 2PN SS 


84.6 


85.6 


83.8 


82.7 



TABLE III: For a typical binary observable by Advanced LIGO, we compare the SNR obtained using the 2.5PN amplitude 
corrected waveform without spin effects, with spin-orbit effects, and with spin-orbit and spin-spin effects. In each column we 
show the component spins (xi -Ln, Xz 'Ln). In all cases we use the SPA phase with non-spinning terms up to 3.5PN order, and 
spin terms up to 2.5PN order, as given in Eqs. (|6.22|) - (|6.25p . The binary is at a distance of 100 Mpc with the same orientation 
as in Table U] 





(1,-1) 


Ui 
(lO'^ 
(0.1,-1) 


sA SNR 
+ 10") A/0 
(0.01,-1) 


(0.001, -1) 


2.5PN 

2.5PN + 1.5PN SO 

2.5PN + 1.5PN SO + 2PN SS 


2532.9 
2909.7 
2924.1 


2532.2 
2545.3 
2546.7 


2525.1 
2503.5 
2503.7 


2576.9 
2551.1 
2551.1 



TABLE IV: For a typical binary observable by LISA, we compare the SNR obtained using the 2.5PN waveform without 
spin effects, with spin-orbit effects, and with spin-orbit and spin-spin effects. In each column we show the component spins 
(Xi ■ Ln, X2 • Ln). In all cases we use the SPA phase with non-spinning terms up to 3.5PN order, and spin terms up to 2.5PN 
order, as given in Eqs. (|6.22|| - (|6.25|l . The binary is at a distance of 3 Gpc with the same orientation as in Table IT] The binary 
masses and distances refer to the redshifted quantities. 

orders remains the same for all the random realizations except between the Newtonian and 0.5PN order. Though on 
most of the occasions, the SNR increases from Newtonian to 0.5PN order, there arc cases when it decreases, albeit 
slightly. All these cases where the SNR decreases have inclination angle 9 very close to zero or tt. For these cases, 
the third harmonic, which is proportional to sin0, is largely suppressed and the spin-dependent interference accounts 
for the small drop in SNR. This drop is observed for systems for which Xi = X2 = — 1 whereas the non-spinning and 
Xi = X2 = 1 cases consistently showed the increase in SNR between Newtonian and 0.5PN order. To further assert 
this, we fix the inclination angle to a value very close to zero and tt and randomly varied the other three angles. We 
find that for all the realizations the SNR decreases in going from Newtonian to 0.5PN order. In brief, the trends 
shown in Table |T] is quite general except for inclination angles close to zero or tt. We however note that the trends of 
Table |T] need not be same for much higher masses when the leading harmonic approaches the lower cut-off frequency 
of the detector (2Flso — /«)■ We have not done a thorough analysis for the whole mass range. 

VII. CONCLUSIONS 

The ongoing search for GWs from compact binaries with the network of interferometers LIGO, Virgo and GEO, 
and the work at the interface between analytical and numerical relativity aimed at providing accurate templates for 
the search, has made it urgent to include higher-order PN effects in the theoretical predictions of the waveforms. This 
paper is a step forward in this direction. 

We provided ready-to-use time-domain waveforms for spinning, precessing binaries moving on nearly circular orbits 
through 1.5PN order and decompose those waveforms in spin- weighted —2 spherical harmonics [see Appendices El 
andlB], Neglecting radiation-reaction effects and assuming ^ L, we found that the inclination angle t between the 
total angular momentum and the Newtonian orbital angular momentum (sec Fig.[T|) is a 0.5PN correction. Motivated 
by this, we expanded the GW polarizations and spin-weighted spherical harmonic modes in a Taylor series in t 
[see Eqs. p.l6p . (|3.17p and Eqs. (|4.16a|) - (|4.16rp ]. Their expressions become much simpler and allow one to extract 
interesting physical features of the gravitational waves from precessing binaries. 

We found that, in contrast to what happens in the non-spinning case, the hi^s are not in general proportional to 
e~™'*. They also depend on the angles t and a, where t is the inclination angle of the Newtonian orbital momentum 
relative to the total angular momentum and a is the angle between the x-axis and the projection of the Newtonian 
orbital angular momentum onto the x-y plane (see Fig. [1]). For example, the terms independent of l are proportional 
to g-""(*+")^ the terms that are linear in l are proportional to e-»("'+™ )(*+«) cos^ or e-»(™+™ )(*+«) sin'J/, while 
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higher-order contributions in l involve terms of the form e~*(™+™ cos" 5" sin^^*, where a,6 £ N and m' € 

— 1,0, 1. In the presence of precession, the angles t and a vary in time and the different harmonics present in each 
of the modes interfere, causing a strong modulation of the mode amplitudes. Wc also found that, in contrast to 
what happens in the non-spinning case, the signal can be largely distributed among modes m) other than the 
(2, 2) mode. With our choice of the source frame, when spins are maximal and the binary system has significant 
mass asymmetry and/or a large inclination angle, we found that the amplitude of the (2, 0) and (2, 1) modes can be 
comparable to the amplitude of the (2, 2) mode, especially during the last stages of inspiral. For the mass ratios we 
considered, wc found that the € = 3 and i — A modes are generally one or two orders of magnitude smaller than the 
£ = 2 modes. These results are summarized in Figs. [2l [3]and|4l for binaries with mass ratio 1:1 and 4:1, and for two 
maximal spin configurations having a small or large inclination angle 6. The ready-to-use time-domain waveforms 
for spinnin g, p recessing binaries can be employed for accurate cornparisons with numerical simulations of binary BHs 
f63l. [79I [qoI. loil. [93l. Im| and for designing time-domain [gl, [gl, 113, HI] analytical templates. 

Restricting ourselves to spinning, non-precessing binaries, we computed ready-to-use frequency-domain waveforms 
in the stationary-phase approximation. We derived IPN and 1.5PN order spin-orbit effects, and 2PN order spin- 
spin (spin(l)-spin(2) only) effects in the frequency-domain GW amplitude [see Eq. (j6.14p . and Eqs. (jDip - ()D2ip in 
Appendix |D]. We also calculated the 2PN spin-spin {including spin(l)-spin(l) and spin(2)-spin(2) effects), and the 
2.5PN order spin-orbit effects in the frequency-domain GW phase [see Eqs. (|6.22p . (|6.25p and (|6.24p ]. For the 2PN 
spin-spin terms, we found agreement with Refs. (3^ . [33| . We wrote the frequency-domain waveforms in a rather 
compact way, so that they can be easily used for data analysis and for building analytical frequency-domain [^, [t^ 
templates. 

In the non-precessing case, we found that, through 2PN order, spin effects in the amplitude affect only the PN 
corrections to the first and second harmonics. Thus, through 2PN order, spin effects do not yet extend the mass 
reach of GW detectors. However, as seen in Figs. [51 [HI they can interfere with other harmonics and, depending 
on the spin orientation, lower or raise the signal-to-noise ratio of ground-based (see Tables HI IIIip and space-based 
detectors (see Tables HIl IIV|) . Wc also expect that those spin terms will help in localizing the binary source in the 
sky. We leave to a future pub lication the use of the waveforms derived in this paper to extend parameter-estimation 
predictions [H, S, Si, [H, 113, [H, [4§, HO, HH, [H, [H, [Hi, i55|, ^5^, H^, [H, [s^ of ground-based and space-based detectors 
to spinning, precessing binaries. 
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APPENDIX A: READY-TO-USE GRAVITATIONAL- WAVE POLARIZATIONS FOR PRECESSING 
BINARIES ON NEARLY CIRCULAR ORBITS THROUGH 1.5PN ORDER: GENERIC INCLINATION 

ANGLES 

In Sec, mil we wrote the GW polarizations 
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expanded in the inclination angle l. Here we give the full expressions. The Newtonian, 0.5PN and IPN order 
coefficients were computed explicitly in Ref. [34| . the 1.5PN order coefficients, are computed for the first time in this 
paper. They read 
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APPENDIX B: GRAVITATIONAL- WAVE MODES FOR PRECESSING BINARIES ON NEARLY 
CIRCULAR ORBITS THROUGH 1.5PN ORDER: GENERIC INCLINATION ANGLES 

In Sec. IIVI we wrote the gravitational-wave modes expanded in the incUnation angle l. Here we give the full 
expression of the /122, ^33 and /121 modes. They read 
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APPENDIX C: CENTER-OF-MASS ENERGY AND GRAVITATIONAL- WAVE ENERGY FLUX 

For nearly circular orbits, the centcr-of-mass energy is known through 2PN order, when spins arepresent and 3PN 
order when spins are neglected. The coefficients entering Eq. (pT^ are [11, [11, [11, III, [8l|, [sl, [11, [13, [111 
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The GW e nerg y flux is known through 2.5PN order for spin effects [23, [30, [S^, [33| , and 3.5PN order when spin effects 
are absent [23|. The coefficients in Eq. (|6.19p read 
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APPENDIX D: FREQUENCY-DOMAIN AMPLITUDE CORRECTIONS 



We give here the complex coefBcicnts Cjl" appearing in the frequency domain non-precessing waveform (|6.13[) . The 
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